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Abstract 

Light quark masses can be determined through lattice simulations in regular- 
ization invariant momentum-subtraction(RI/MOM) schemes. Subsequently, 
matching factors, computed in continuum perturbation theory, are used in or- 
der to convert these quark masses from a RI/MOM scheme to the MS scheme. 
We calculate the two-loop corrections in quantum chromodynamics(QCD) to 
these matching factors as well as the three-loop mass anomalous dimensions 
for the RI/SMOM and RI/SMOM^ schemes. These two schemes are char- 
acterized by a symmetric subtraction point. Providing the conversion factors 
in the two different schemes allows for a better understanding of the system- 
atic uncertainties. The two-loop expansion coefficients of the matching factors 
for both schemes turn out to be small compared to the traditional RI/MOM 
schemes. For rif = 3 quark flavors they are about 0.6%-0.7% and 2%, respec- 
tively, of the leading order result at scales of about 2 GeV. Therefore, they will 
allow for a significant reduction of the systematic uncertainty of light quark 
mass determinations obtained through this approach. The determination of 
these matching factors requires the computation of amputated Green's func- 
tions with the insertions of quark bilinear operators. As a by-product of our 
calculation we also provide the corresponding results for the tensor operator. 



1 Introduction 



Light quark masses, like the up-, down-, and strange-quark masses, are fundamental 
parameters of the Standard Model, and their precise determination is thus an impor- 
tant task. They can be determined, for example, with the help of lattice simulations in 
combination with nonperturbative renormalization(NPR). In this context regularization 
invariant momentum-subtraction (RI/MOM) schemes [1] play a crucial role; for a recent 
overview see e.g. Ref. [2]. 

In Ref. [3] these light quark masses were determined in the RI/MOM scheme and 
subsequently converted to the MS scheme (HE]. This conversion requires the computation 
of a matching factor C^, which performs this transformation of the quark mass from 
the scheme x into the MS mass. Since the RI/MOM schemes do not depend on the 
particular regulator which has been used to regularize the ultraviolet divergences, this 
matching factor can be calculated in continuum perturbation theory using dimensional 
regularization [B]. In the RI/MOM scheme this matching factor is known up to three-loop 
order in perturbative QCD [UlTllH]; the same holds for the RI'/MOM [8l[9] scheme. 

Both the RI/MOM and RI'/MOM schemes employ a renormalization procedure which 
uses an exceptional subtraction point. Here the subtraction point defines the configuration 
of the external momenta of the considered amplitude, where the ultraviolet divergences 
are subtracted. In continuum perturbation theory, mass and fermion field renormalization 
constants are typically computed by considering higher order corrections to fermion self- 
energy diagrams. The fermion self-energies are related through Ward-Takahashi identities 
to amputated Green's functions of quark bilinear operators, i. e. the vector, axial- vector, 
scalar and pseudoscalar operators. These relations allow an extraction of the renormal- 
ization constants from these amputated Green's functions with operator insertions, rather 
than from self-energies. In the case of an exceptional momentum configuration, no mo- 
mentum transfer leaves the operator. 

However, at low scales (~ 2 GeV) the perturbative expansion of the matching factors 
for the RI/MOM and RI'/MOM schemes exhibit a poor convergence behavior and, as a 
result, introduce large systematic uncertainties in the determination of the MS masses 
of the light quark sector. This amounts to approximately 60% of the total error. Apart 
from that, the lattice simulations are, for an exceptional subtraction point, more prone to 
unwanted infrared effects. For this reason, the use of a symmetric subtraction point was 
proposed in Ref. [10], and the concepts and framework of these new RI/SMOM schemes 
have been worked out in Ref. [TT]. A nonexceptional or symmetric subtraction point is 
characterized by the fact that a momentum leaves the operator of the amputated Green's 
function, and in the case of a symmetric subtraction point, the squares of all momenta 
leaving the amplitude are equal. A nonperturbative test of such a RI/SMOM scheme can 
be found in Ref. [12]. The one- loop QCD corrections of these matching factors in the 
RI/SMOM schemes are known [TT], and are shown to have a better convergence behavior 
than the traditional RI/MOM and RI'/MOM schemes. However, the question of whether 
this behavior persists at higher orders in perturbation theory still remained unanswered. If 
it were confirmed, this would lead to a significant reduction of the systematic uncertainties 
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associated with the determination of the hght quark masses. 

The goal of this paper is to extend the work of Ref. [11] and to provide the two-loop 
QCD corrections to the matching factors of the so-called RI/SMOM and RI/SMOM^^ 
schemes as well as the three-loop mass anomalous dimensions. 

This paper is structured as follows: In Section [2] we introduce the notations and con- 
ventions that are used. In Section [3] we give an outline of the perturbative calculation, 
and in Section H] we present the results for the matching factors and the mass anoma- 
lous dimensions. Finally, in Section [5] we close with a summary and conclusions. In 
the appendixes we provide additional information about some master integrals and, for 
completeness, results for conversion factors with complete gauge parameter dependence 
as well as results for anomalous dimensions. 



2 Generalities and Notation 

In Ref. [H] two regularization invariant momentum-subtraction schemes with a sym- 
metric subtraction point have been defined, the so-called RI/SMOM and RI/SMOM^^ 
schemes. The RI/SMOM (RI/SMOM^^) scheme can be seen as an extension of the 
RI'/MOM(RI/MOM) scheme from an exceptional to a nonexceptional subtraction point. 
The renormalization conditions of these two new schemes are given in Eqs. (10), (11), and 
(A15) of Ref. [11] and allow for a determination of the renormalization constants of the 
fermion field and the quark mass m through the computation of the nonsinglet ampu- 
tated Green's function Aq with the insertion of quark bilinear operators O, that is to say, 
the scalar {O = S), the pseudoscalar (O = P), the vector {O = V), and the axial- vector 
(O = A) operators. Further details can be found in Ref. [TT] , whose conventions we follow 
for these quantities. 

In the determination of the light quark mass through lattice simulations, the RI/SMOM 
and RI/SMOM^^ schemes serve as intermediate schemes before the conversion of the quark 

mass to the MS scheme. This conversion is performed in continuum perturbation theory 
through the computation of a matching factor with the property 

= C^m^, X G {ri/smom, Ri/SMOM^^}. (1) 

The index r denotes here and in the following a renormalized quantity. As shown in 
Ref. [H] these conversion factors can be extracted from the amputated Green's function 
of the pseudoscalar (or scalar) operator through the equation 



X G {rI/SMOM, RI/SMOM^ }, (2) 



7m ^ 

sym 



where is the corresponding matching factor for the fermion field which converts the 
field from the RI/SMOM(RI/SMOM^^) scheme to the MS scheme via the equation 

= VQ^I' X G {ri/smom,ri/smom^^}. (3) 
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The matching factor can be determined through the computation of the amputated 
Green's function with the insertion of the axial- vector (or vector) operator 



(4) 

sym 



and 



(5) 



RI/SMOM 1 .;.,MS. ^ 

"^^o 48 L ' J 
where Eqs. and (E]) are distinguished through the use of two different projectors in the 
trace with which the amputated Green's functions are muhiphed. The subscript sym in 
Eqs. ([2]), (jl]) and ([5]) stands for the restriction of the amputated Green's function to the 
symmetric momentum configuration 

f{pl,pl,q'^)\,y^ = f{pl,pl,q'^)\p2^p2^g2^_^2 , (6) 

where pi and p2 are the momenta of the external fermions and q = Pi —p2 is the momentum 
transfer leaving the operator. The momenta are also defined in the two-loop diagrams of 
Fig. [H The symbol /i^ > denotes the renormalization scale, which we choose for the 
RI/SMOM(RI/SMOM^^) scheme to be equal to the one in the MS scheme. 

Also, the tensor operator {O = T) has interesting applications in the context of lattice 
simulations; see e.g. Ref. [13]. Following the notation and Eq. (12) of Ref. [11] one 
obtains the corresponding matching factor C^, which converts the tensor operator from 
the X scheme to the MS scheme, 

OT = C^oOw (7) 

by computing 



C^ = C^ lim -i-Tr 

^ 9 mR^O 144 



X e {rI/SMOM, RI/SMOM }, (8) 



7m- 

sym 



with cr^^ = |[7/.,7;.]- 

The conversion factors in Eqs. ([T]), (E]), and are gauge dependent. In the case 
of the mass conversion factor this gauge dependence compensates for the one in the 
mass renormalization constant which is obtained through lattice simulations; the NPR 
procedure of Ref. [1] is gauge dependent. The lattice simulations are typically performed 
in the Landau gauge (^ = 0). For generality we keep the complete dependence on the 
gauge parameter and use for the gluon propagator 

^(-«-H-(l-04^). (9) 

3 Calculation 

The calculation of the two-loop QCD corrections proceeds in two steps. In the first step, 
after the generation of the diagrams, all loop integrals that appear are mapped on a small 
set of master integrals. In the second step these master integrals need to be solved. 
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Figure 1: Two- loop Feynman diagrams which contribute to the computation of the non- 
singlet, amputated Green's functions. The spiral lines denote gluons, the solid lines 
represent fermions, the dashed lines are ghost fields, and the black boxes indicate the 
inserted operator. 



3.1 Manipulation of loop integrals 

In order to generate the required Feynman diagrams, shown in Fig. [H we have used the 
program QGRAF [T3]. We identify the different topologies and adopt the proper notation 
with the help of the packages q2e and exp [TSHTT] . This allows us to prepare an output 
which can straightforwardly be used to perform the reduction to master integrals. The 
reduction is achieved with the traditional integration-by-parts (IBP) method [18] in com- 
bination with Laporta's algorithm [T^I2U] and has been carried out with a FORM pTH25] 
based implementation. The rational functions in the space-time dimension d, which ap- 
pear while solving the arising linear system of equations, have been simplified with the 
program FERMAT pi] . 

In the case of the computation of the amputated Green's function with the insertion 
of the pseudoscalar or axial-vector operator, we use a naive anticommuting definition of 
75 for the treatment of 75 in dimensional regularization [SlESj, which is a self-consistent 
prescription for the flavor nonsinglet contributions considered here [26||27]. 



3.2 Master integrals 

After the IBP reduction of the two-loop amplitude, seven massless master integrals sur- 
vive. The one- and two-loop master integrals are shown in Fig. |2l They are defined in 
Minkowskian space in d = A — 2e space-time dimensions by 
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Figure 2: There appear two one-loop master integrals and seven two-loop master integrals 
from which two are factorized. The symbol A^pr denotes a /-loop topology with p lines. 
The number r serves as a running number to enumerate topologies with the same number 
of loops and propagators. 
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with the denominators 

Di = ij + le, D2 = {ii + gs)' + ^e, D, = {h - q^f + ze, 

D^ = ll + ze, = (£2 - q2f + ^e, = (£1 - £2)' + ^e, 

^7 = (4 - ^2 - gi)' + ^e, 

where li and £2 are loop momenta and gi, 52? and are external momenta with q\ = 
I2 — I3 — ~/^^- The symbol e ~ 2.71828 is Euler's number and 7^ ~ 0.577216 is the 
Euler-Mascheroni constant. 

The one- and two-loop self-energy integrals A^2i''' -M3? are well known to all orders 

f2) (2) 

in e. The factorized master integrals and A4l,i can be obtained from taking the 

(2) (2) 

product of the two one-loop integrals. Each of the master integrals M.42 and can 
be written as a product of a scalar one-loop two-point function and a one-loop three- 
point function A^3i^ , which has a noninteger power of one of the propagators. The latter 
has been determined analytically in terms of hypergeometric functions to all orders in 
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e in Ref. [28]. Algorithms to expand generalized hypergeometric functions in a small 
parameter, like e, have been developed in Refs. [29l[30]. For the practical implementation 
we use the Mathematica packages HypExp and HPL [3TI433] . Expansions of hypergeometric 
functions can also be found in Refs. [MlEEj. The master integrals M.\,2 A^gi 
be taken from Ref. [SB] in the special case of symmetric external momenta = — /^^ for 
i = 1,2, 3. Results for the master integrals in terms of harmonic polylogarithms can also 
be found in Ref. [37] . 

The results expanded in e to the maximally required order read 
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where \E'(x) is the digamma function \E'(x) = r'(x)/r(x); here the prime denotes the 
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derivative of the T function. We also define Sn{x) = -^Im 



with the poly loga- 
rithm function Lin{z) = Yl'^^i f^- The symbol C„ = Li„(l) is the Riemann zeta function. 
The constants T-L^i and T-L^^ can be expressed in terms of harmonic polylogarithms and 
are discussed in Appendix El 

The higher orders in e of some of the master integrals are needed due to the appearance 
of so-called spurious poles which arise while solving the linear system of IBP equations. 
An alternative approach is to find a different set of master integrals through the method 
of e-finite basis [38] which exploits the fact that the choice of the master integrals is not 
unique and which requires only the evaluation of the master integrals up to the order e". 

The master integrals constitute an essential input of our calculation, so we have 
checked the results using traditional Feynman parametrizations, partially analytically and 
partially through the numerical evaluation of the resulting Feynman integrals in order to 
ensure their correctness. 



4 Results 

4.1 Matching factors 

In this section we give the conversion factors up to two-loop order in perturbative QCD for 
the fermion field and the mass parameter. This is accomplished as previously described in 
Section |2l For this purpose we decompose the matching factors into a part which depends 
on the gauge parameter and a part which is gauge parameter free, 

= ^y,L + C*^,?' ^ ^ {ri/smom, ri/smom^^}, y € {q, m}, (19) 

where Cy^ contains all the terms which depend on the gauge parameter ^, whereas ^ is 
the conversion factor in the Landau gauge. In the following we present the contributions 
in the Landau gauge; for completeness, we also provide the gauge dependent terms in 
Appendix [Bl The one-loop results have been determined in Ref . [11] . 

Let us start with the conversion factors for the fermion field of Eqs. and ([5]). They 
are important for the light quark mass renormalization procedure in Eq. ([2]), but also 
enter in the renormalization of other interesting quantities, i.e. any multiquark operator 
in such schemes, like, for example, the parameter. With the help of Ward-Takahashi 
identities it has been shown in Ref. [TTI that (^^Vsmom _ ^ri /mom^ have checked 



that our order result for (^^Vsmom agreement with C^^ /mom ^^^^ j^^^ jg ^ rpj^^ 
fermion field mate 
result is given by 



fermion field matching factor in the RI/SMOM^ scheme at two-loop order is new. The 
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The symbols Cp and denote the Casimir operators of the fundamental and adjoint 
representations of SU(A''c). The normalization of the trace in the fundamental represen- 
tation of SU(A'^r;) is given by Tp = 1/2. For the number of colors A'^c = 3 one obtains 
Cp = 4/3 and Ca = 3. The number of active fermions is given by n^. 

Starting from two-loop order the strong coupling constant also has to be renormal- 
ized, and it has to be ensured that the Slavnov- Taylor identities, which provide relations 
among the different renormalization constants, are preserved. These relations guarantee 
that the renormalized coupling constant remains the same for the different interaction ver- 
tices of the QCD Lagrangian. Therefore, we perform the renormalization of in the MS 
scheme with a^,/? = Z~'^^~'^'^as^B- When not considering the Landau gauge, we also have 
to perform a renormalization of the gauge parameter = Z7^^b with = Z3, where 



a,fj, 



the renormalization constant Z3 is the one of the gluon field C^'^ with G^'^ = 1/ y/ZsG^ . 
We renormalize these quantities in the MS scheme (see also Appendix [B]) . In general, the 



ghost field c"" also needs to be renormalized, = I/1 



Z^c°Q^ which can also be done in 



the MS scheme, but is not needed within this calculation. 

In complete analogy to Eq. fl20|) . the mass conversion factor, required for the matching 
of the light quark masses determined in the RI/SMOM scheme to the MS scheme, reads 
up to two-loop order in the Landau gauge = 0) 
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where the symbol S is given in Appendix |Al Similarly, we obtain, for the mass conversion 
factor up to two-loop order in the RI/SMOM„, scheme, 
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We want to mention that the matching factors in Eqs. (12T|) and (122!) are related to con- 
version factors of quark bilinear operators. As shown in Ref. [H] Ward-Takahashi identi- 
ties allow us to write the matching factor Cf,, which converts the pseudoscalar operator 
renormalized in the x scheme to the MS scheme, as the inverse of for x G {ri/smom, 
Ri/SMOM^^}. Furthermore, the latter are also related to the corresponding matching factor 
of the scalar operator C'§ = Cf, = 1/C^^, x E {ri/smom, ri/smom^^}. For the definition 
of these matching factors we adopt the same notations as in Ref. [11]. For completeness 
we give here also the result for the matching factor of the tensor operator and again 
disentangle into its Landau gauge (^ = 0) and its gauge dependent parts in complete 
analogy to the decomposition of Eq. ( IT9|) for y = T. The Landau gauge component is 
given by 

~ . ,,(-)c,(i*'(i)J,^-l- 



a 



9 



/a,\2 r [473 20 ^,.1. 40 2 



599 l!!^.m_l^,m._ 356 16 n 
9 ^ 27 27 ^ 81 



24 



1^ TT- 



32 
243' 



-TT 



--vE^'"(l) - 24052(f) + 120^2(1 ) + 16053(f) - 20053(1) 



-CaCf 
4 



10 log (3) TT 51og2(3)7r MSvr^ 



8491 755 
" 216 ~ 108 

243^ 24 

16 log (3) 7r 



6^3 
2 

27 



162v^ 



I2C3 



755^ 
162' 



-TT 



1285,(1)- 6452(f) -^53(f) 



3^3 



41og2(3)7r 1167r3 50 
9v^ ^ 24373 ^ 



81 

256 . , 320 
— ^3(f) + — ^3(1) 

+ (23) 



0.21517295 



47r 



^ j [43.38395007 - 4.10327859^/] + 0{al 



and 

RI/SMOM. 



c, 



-l/J. 



1 + 
+ 



+cl 



Air 

Att) 
183 



3 



*'(!) 



2 . 1 



9 3 
299 _ 20 
54 27 ^''> 



40 



-TT 



245 

H ^ (^) 

18 



27 



245 



-TT 



TT 



38 
243' 



-TT 



^'"(1) - 19252(f) + 9652(f) + 12853(f) - 16053(f) + S 



+CaCf 



108 
81og(3)7r 

3185 



21og2(3)7r 587r3 



3^3 



108 



81^3 
n2 



14C3 
193 



14 



-TT 



81 243 



1 



49 log (3) TT 491og2(3)7r 142l7r=^ 



38 



12^3 



144^3 3888^3 3 



C3 



(24) 



Atx 



1.11816038 



47r 



[8.607630493 - 1.955130440n/] + C»(af). 



The gauge dependent component can be found in Appendix [Bl The one-loop order in the 
RI/SMOM scheme has already been computed in Ref. [11]. The result in the RI'/MOM 
scheme is available up to three-loop order in Ref. [9]. 
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4.2 Comparison with results using exceptional momenta 

In order to analyze the quality of the new two-loop order of the mass conversion factors 
in Eqs. (12T]) and (l22l) . we compare them to the results of the traditional RI'/MOM and 
RI/MOM schemes. For this purpose we evaluate the matching factors numerically, again 
in the Landau gauge, for Uf = 3 flavors. For the strong coupling constant we insert 
ajn = 0.1 0. The results in the RI/SMOM versus the RI'/MOM scheme for A^c = 3 read 

^Ri/SMOM ^ 1 - 0.0161380... -0.00660442..., (25) 
^RiVMOM ^ 1 - 0.1333333... -0.07585848... - 0.0556959..., (26) 

where each term stands for the next order in perturbative QCD, i.e. the leading, next- 
to-leading, and next-to-next-to-leading order results. Similarly follow the results for the 
RI/SMOM^^ versus the RI/MOM scheme 

RI/SMOM^ , ^ 

C^'l = 1 - 0.0494713... - 0.0228421..., (27) 

^RI/MOM ^ 1 - 0.1333333... - 0.0815876... - 0.0602759.... (28) 

Both the RI/MOM and RI'/MOM schemes are known to three- loop order, where we 
have taken the results of Ref. [8]. One can observe that the matching factors of the 
schemes with a symmetric subtraction point in Eqs. (1251) and (1271) show a much better 
convergence behavior than the ones in the RI'/MOM and RI/MOM schemes of Eqs. fl26l) 
and fl28l) which are characterized by an exceptional momentum-subtraction point. This 
observation has already been made in Ref. [H] at one-loop order and is now confirmed 
at the two-loop level. The size of the three-loop corrections is about 6% of the leading 
order result for both the RI/MOM and RI'/MOM schemes, whereas in the RI/SMOM 
scheme the two-loop result is already significantly smaller, being of about only 7 per mili 
of the leading order one. Similarly, for the RI/SMOM^ scheme the size of the two-loop 
term is of about 2%, which is also smaller than the three-loop terms of the two schemes 
with an exceptional subtraction point. The use of these new results in light up-, down-, 
and strange-quark mass determinations in the context of NPR will allow us to reduce 
the uncertainties related to the matching procedure to the MS scheme, due to a smaller 
truncation error of the perturbative series, which will significantly reduce the error of 
these light quark masses obtained in this approach. 



^Thc choice of this value corresponds approximately to the value of a^/Tr at a scale of 2 GeV and is 
sufficient to illustrate the effect of the new results. It allows the reader to easily replace it by her/his own 
value of as. 

^ Since the percentage correction of the two- loop term in the RI/SMOM scheme is very small, its size 
is more sensitive to the exact value of than in the other schemes. Thus a slightly smaller value of 
tts/TT leads to a slightly smaller correction of about 0.6% for the RI/SMOM scheme. 
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4.3 Anomalous dimensions 



The mass anomalous dimension of the RI/SMOM and RI/SMOM^^ schemes is used to 
run the quark masses to different energy scales. It is defined by 



7„ 



dlog rrf{ii) 



-7n 



Aw 



Aw 



An/ 



0{at 



(29) 



d log /i^ 

where the superscript x stands, similar to Eq. (fT9l) . for the scheme with x G {ri/smom, 
Ri/SMOM^^, Ms}. Both mass anomalous dimensions are related to the MS mass anomalous 

dimension through the conversion factors C^^ presented in Section 1411 (see Ref. [8]). 
In the Landau gauge holds 

X e |ri/smom, ri/smom^^I, 



' m 



-(3- 



(30) 



where the QCD /3 function and are given, for completeness, in Appendix O The 
one- loop mass anomalous dimensions in the RI/SMOM and RI/SMOM^ schemes are 

equal to the one in the MS scheme; the mass anomalous dimension at two-loop order has 
been determined in Ref. [11] and is also given in Appendix [Dl The three-loop result is 
presented here and reads, in the Landau gauge in the RI/SMOM scheme|l 



_ (2),RI/SMOM^ 
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--vp'd) + 88.2(f) - 176.2(f) - —.3(1) + —33(f) + — C3 



+nfTF { C 



77 



76 
320 



1527r2 

27 
2727r2 



1167r3 
243^3 



4647r^ 
243 



8 y/yn I 47rlog^(3) 

27^^ 



f ^'(i) + ^^'(i) - ^'i^ - + 64..(|) - 128..(f ) 



.,.256 176 
^^3(3) + — ^3(2)-— C3 



3 



87rlog(3) 27rlog2(3) 



+ CaCf 
1534 



3^3 
7870 



27 



30687r2 587r3 
—r-. — + 



567r^ 



243^3 243 



647r2_^.,,,, 16..,., ,2 



3^3 



9^3 



+ _ 32.2(f) + 64.2(f) + ^.3(f ) - ^.3(f ) + ^C3 



27 
L6C 

y 



^'(1) -i^^'d) + §*'(!)' 



+CaC\ 



-9 



6167r2 
27 



3197r3 
243^3 



12767r^ 
243 



22s + 



3 ^^'^ 3 
447rlog(3) ll7rlog^(3) 308 



3^3 



9^3 



9 



■^The minus sign on the left-hand side of Eqs. ((3T|) . (|32p and (|35|) is correct, but was missing in the 
first arXiv version. It is needed to be consistent with the definition in Eqs. (f29|) and (p3|) . respectively; 
see also the note added. 
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H- - - 17fa,(|) H- 352,,(|) 

856 3207r2 160 



^80 



^3(1) 



704 



«3(f)+^C3 



27 



27 



In the RI/SMOM scheme the two-loop order is again available in Ref. [TT] (see also 



Appendix 0) and the three-loop result in the Landau gauge is given by^ 
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6 W 



121 
~9~ 



81 



v|/'(l)2 + 44^3 + + 22S - 

27 ^ 243 54 Vs; 



67715 42867r2 
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243 
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^3(1) 

10 
27 



2143 SSvr^ _ 22 
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-Cs 



80 
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'27 



^'(l)'-64C3 



+CaCf 
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^ 81 
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27 



(32) 



Similarly to the mass anomalous dimension we define the anomalous dimension of the 
fermion field as 



la 



d log yU^ 



(0),x^ 



7, 



(2),x 

in J 



0{at 



(33) 



where the renormalized and bare fields are connected by = -^Z^i^. The relation 
between the MS fermion field anomalous dimension and the ones found in the RI/SMOM 
and RI/SMOM^^ schemes is, in the Landau gauge, again given by [8] 



Iq = - f^-^ir^ ^ ^ {ri/smom, ri/smom^^,ms}. 



(34) 
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In the RI/SMOM scheme ^^^/^'^^'^ jg equal to /mom loop orders in perturbation 
theory, which has been shown in Ref. [11]. We have reaffirmed this equivalence to 0{a^) 
by comparing our results with those of Ref. [8]. For the sake of brevity, we omit this 
result since 7^^ /mom^ three-loop level, can be found in Ref. [8]. However, such an 

equality is no longer true in the RI/SMOM^^ scheme. Therefore, we provide the three- 
loop contribution to the fermion field anomalous dimension in this scheme in the Landau 
gauge below^ 
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160 
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2^ ^280 
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(35) 



The one- and two-loop contributions are again given in Ref. [TT] and can also be found in 
Appendix [Dl 

The anomalous dimension of the tensor operator is defined by 

2 /„\ / rv . \ 3 



dlogZ^ 

It ^ 



d log /i^ 



(0),x «s {l),x ( OLs 



It 



47r 



(36) 



where the renormalized and bare operators are related by the renormalization constant 
with 0\ = Ob ioT O = T and x G {ri/smom, ri/smom^ }. The anomalous dimensions 
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can again be obtained to three-loop order in the Landau gauge by inserting the result for 
the matching factors 2, of Eqs. (123|) and (124|) into 



It = 



(3 



d log C; 



T,L 



X G {ri/smom, ri/smom^^} 



(37) 



The MS anomalous dimension 7^?^ of the tensor operator is known up to three-loop order 
and can be taken from Refs. |39lll0l|9]. It is also given explicitly in Appendix O The 
results for the RI/SMOM and RI/SMOM^^ schemes read 
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5 Summary and Conclusion 



We have computed the two-loop QCD corrections for matching factors which convert light 
quark masses renormalized in the RI/SMOM and RI/SMOM^^ schemes to the MS scheme. 
These schemes are extensions of the traditional regularization independent momentum- 
subtraction schemes, like the RI/MOM or RI'/MOM scheme. The latter are characterized 
by an exceptional subtraction point, whereas in the RI/SMOM and RI/SMOM schemes 
the ultraviolet divergences are subtracted at a symmetric subtraction point. This allows 
for a lattice simulation with suppressed contamination from infrared effects. The per- 
turbative expansion coefficients at two-loop order are significantly smaller than in the 
traditional RI/MOM and RI'/MOM schemes, with about 0.6%-0.7% and 2% of the lead- 
ing order result for the RI/SMOM and RI/SMOM,^ schemes for Uf = 3 flavors at scales 
of about 2 GeV. We also provide the results for the three-loop anomalous dimensions of 
the quark field and the mass in both schemes. The mass anomalous dimension can be 
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used to run the quark masses to different energy scales. 

The use of these matching factors will reduce the uncertainties associated with the 
matching procedure for converting the quark mass from the regularization invariant 
momentum-subtraction scheme to the MS scheme due to smaller truncation errors of the 
perturbative series. These results will allow for an MS light quark mass determination 
only together with lattice simulations in the context of nonperturbative renormalization, 
with a significant reduced systematic error compared to previous determinations which use 
momentum-subtraction schemes with an exceptional subtraction point, whose matching 
factors show a poor convergence behavior at the required energy scales. 

In addition, the availability of both the RI/SMOM and RI/SMOM^ schemes at two- 
loop order will also allow for a better assessment and control of the uncertainties. This 
can be achieved by converting results derived in the different regularization invariant 
momentum-subtraction schemes to the MS scheme, using the above matching factors, 
and subsequently comparing and cross-checking the obtained MS results. 

Our computation is accomplished by studying at two-loop order quark bilinear oper- 
ators inserted into amputated Green's functions for the vector, axial-vector, scalar, and 
pseudoscalar operators. We also provide the corresponding results for the matching factors 
and the anomalous dimensions of the tensor operator. 
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results for Cg and Cm the corresponding anomalous dimensions agree 

up to a global sign which has been rectified in the current version. We would like to thank 
the authors of Ref. [H] for communications and for pointing this out to us. 



A Master integrals and harmonic polylogarithms 

The harmonic polylogarithms H of Ref. [12] are defined recursively as 

H{a,ai,. . . ,ak;x) = dx'fa{x')H{ai,...,ak;x') (40) 
Jo 

with 

fi{x) = i^, Mx) = i, f-i{x) = 
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and 

= -log(l -x) , H{0]x) = log{x) , =log(l + x). 

The constants and 7/43 in the master integrals of Eqs. (fTTj) and (fT6|) can be expressed 
in terms of harmonic polylogarithm^; their numerical evaluation is given b}|§ 

= -6.11477558..., = +12.45994893.... (41) 

In the conversion factors these constants always arise in the form of the sum S = T-L^-{ + 

(2) 

1-L\^ in which some of the harmonic polylogarithms canceL It can be expressed with 
^0 = ^ by 

S = -1= [l27rC3-127ri/(0,+,+;^o)-67ri/(-,+,+;2o)+^7r2i7(-,+;^o) 

+2iTX^H{Q, +; Zq) - 18iH{-, +, +, +; ^o) - 36zi/(0, +, +, +; Zq)] (42) 
~ 6.34517334592058543866471875422848044360, 

where the functions H{±, ■ ■ ■ , ±; a;) are defined by the following linear combinations: 

H{+-x) = H{l-x) + H{-l-x), (43) 
H{-;x) = H{l;x)-H{-l;x), (44) 
H{±,ai,...,ak;x) = iJ(l, Oi, . . . , a^; x) ± if (-1, ai, . . . , a^; x). (45) 



B Gauge dependent parts of the conversion factors 

In this section we provide the gauge dependent components of the conversion factors, 
where the gauge parameter ^ has been renormalized in the MS scheme. As mentioned in 
Section 14.11 the fermion field conversion factors in the RI/SMOM scheme are given by 
the ones in the RI'/MOM scheme; therefore we omit them here. However, we did check 
that the gauge dependent terms were also in agreement with the previous calculations in 
the RI'/MOM scheme. The results read as follows: 



RI/SMOM^^ 
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3^3; 2 
e - 1 + 



27r2 



2C3 



3l7r2 



297r3 



47r4 



+24.2(1) -48s2(f) -40s3(f) +32s3(f) 



162^3 27 

27rlog(3) 



*The £ expansion of the master integral m!'-^^ can also be obtained with Refs. j43ll44] . 
^Wc have checked that our numerical evaluation using traditional Feynman parametrizations of these 
constants agrees with the numerically evaluated analytical result. 
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C The QCD /? function and anomalous dimensions 
in the MS scheme 

The QCD /3 function is defined by 



/3 



d log /i^ 

and known up to four-loop order 



. The lowest two orders are given by 



/9i = ^Cl-iCpTFUf -'^CaTfUj. 



(51) 

(52) 
(53) 



The MS mass anomalous dimension [55l|56l[27l[571|58] is defined in Eq. ( I29ll and reads up 
to order 



7: 



(0),MS 



3Cf 



(54) 
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7ii)'^' = \cl+^^C,C^-^^C,T,n,, (55) 

7i^)'^ = - l^C^C^ + Hll^C^Cl + C^T^n, (48C3 - 46) 

/ 55fi\ 14n 



The lowest orders of the MS anomalous dimension of the fermion field [8l[9] defined in 
Eq. fl5^ for ,^ = (Landau gauge) are given by 



if^-""' = 0, ^WMS^_c^ + c,CAf^-2C,T,n,, (57) 
7^ = Cf^- CaCf ^— - 12C3 J + C^Cf - -C3 

+nfTF (SC^ - C^^C^F^) + CFujT'^^. (58) 

The Casimir operators in Eqs. fl52l)- fl58l) of the SU(A''c) group are given by Ca = Nc, 
Cf = {Nl — l)/(2iVc), where Nc is the number of colors. 

For our calculation of the tensor anomalous dimensions in the RI/SMOM and RI/SMOM^^ 
schemes, we need the corresponding MS anomalous dimension of the tensor operator up 
to three-loop order [39ll40ll9]. which we have taken from these references. The first three 
orders read 

(0),MS ^ (l),MS ^2 19 , n n '^^^ n ^c;n^ 

7^^' = Cf, 7t = -Cf^^CfCa^ -CFUfTF — , (59) 



(2), MS 

7t 



F 



9 / - ^ 1 27 



-CFn)Tl. (60) 



D The one- and two-loop anomalous dimensions in 
the RI/SMOM^^ ) scheme 

The one- and two- loop anomalous dimensions have already been determined in Ref. [TT] . 
In order to avoid conventional subtleties we also provide these results here for complete- 
ness. The one-loop results in the Landau gauge are given by 

^,o,,m/SMOM ^ 3C. = 7i°'"""'™^-. ,f = 0, (61) 

^m.RI/SMOM ^ C'^_^l,»)'M/SMOM„^ ^^^^ 



The two-loop results in the Landau gauge read 

3 
'2 



^{i),Ri/SMOM ^ _cf CaCf^^^^ 



9 3 ) 
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(i),Ri/SMOM^^ 3 /69 22 2 

- C^T^Ufl^U'il)- 10 -^n^y (64) 

(l),RI/SMOM 3 ^ ^ 31 „ _ 2 

-7? " = C^-CACF— + CFTFnf-, (65) 



_ (i),Ri/sMOM ^ c'^^l-CAcJ—--^'a) + -A 

'T F 2 ^ ^ \^ 6 9 27 y 



(66) 



(67) 
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